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pi ' Abstract 

, We demonstrate, using the symbolic method together with p-adic and resultant 

methods, the existence of systems with exactly one or two generalized symmetries. 
, Since the existence of one or two symmetries is often taken as a sure sign (or as the 

definition) of integrability, that is, the existence of symmetries on infinitely many 
orders, this shows that such practice is devoid of any mathematical foundation. Ex- 
' tensive computations show that systems with one symmetry are rather common, and 

with two symmetries are fairly rare, at least within the class we have been considering 
in this paper. 

^ 1 Introduction 

> 

CN ■ In 1980 an observation was made at least twice by different authors. In P] it is written 

. Another interesting fact regarding the symmetry structure of evolution equa- 

^) \ tions is that in all known cases the existence of one generalized symmetry 

CN ' implies the existence of infinitely many. 



and in |10| the same statement is made together with the footnote 

^ ' This is not true for systems of equations. For example, the system ut = U2 + 

'i> • v'^/2,vt = 2v2 has a nontrivial group, but this group is exhausted by the one- 

^ ', parameter (with parameter r) group of transformations: Ur = U3 + 3vvi,Vr = 

H : 4i;3. 

Here V2 stands for we use the same notation in this paper. However, the 'counter- 
example' given there is an integrable system, cf 0, section 3. In spite of this fact Fokas 
adapted the remark and formulated the following conjecture in 1987, [Q. 

Conjecture 1 (Fokas). If a scalar equation possesses at least one time-independent non- 
Lie point symmetry, then it possesses infinitely many. Similarly for n-component equations 
one needs n symmetries. 
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Besides their mathematical interest, the observation and conjecture are of some prac- 
tical importance since they are used to argue that it is enough to find only one or two 
symmetries of a system in order to declare it integrable, cf [13| and the discussion in 



10]. This would be reasonable practice if it was simply not possible to prove integrability, 
but the methods employed in [15, ^, |l^, 14, ^ showhow one can effectively obtain inte- 



grability proofs. Moreover, four years later Bakirov |]1[| published the first example of a 
non-integrable equation in the possession of a generalized symmetry. The system 

ut = 5ii4 + vl 

Vt = f 4 

has a sixth order symmetry 

Ut = llue + 5woW2 + 

Vt = Vq 

as one can easily check. It was shown (with extensive computer algebra computations) 
that there are no other symmetries up to order 53. The authors of |Q proved in 1998 that 
the system of Bakirov does not possess another symmetry at any higher order, thereby 
proving that indeed one symmetry does not imply integrability. In [18| it is proved that 



there are in fact infinitely many fourth order systems with finitely many symmetries. The 
method used there could be used for seventh order systems as well and a system with two 
symmetries was found, a counterexample to Fokas' conjecture. However for systems whose 
order is more than seven the approach can no longer be used. The method introduced here 
(using resultants) makes it possible to explore the symmetries of higher order systems. 



2 The symmetry condition 

We study symmetries of 

aiUn + K{v,Vi, . . .,Vn-l) 



Ut 
Vt 



a2Vr, 



where ai,02 € C and K is polynomial in v{x,t) and its derivatives v^- We call these 
systems, where n is the order of the system. 

The right hand side of the equation can be interpreted as an element in a Lie algebra, 
the Lie bracket is computed using Frechet derivatives as follows, see |12], [11| 



aiUn + K 

a2Vn 



biUm + S 
b2Vm 



biD"^ Ds 

>2i 



aiUn + K 

a2Vn 



aaZ?" 



biUm + S 
b2Vm 



aiD^'S - a2DsVn - biD"^K + b2DKU, 




(2.1) 



We call S a symmetry of K if [K, S] vanishes. K is called integrable when there exist 
symmetries on infinitely many orders and almost integrable when there exist symmetries 
on finitely many orders, cf. ll^ 
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How to solve the equation [K, S] = 0, given the order of K and SI First of all, we see 
that if K is polynomial, S has to be polynomial too, cf. [||. This enables us to use the 
symbolic calculus as developed in Q. With the symbolic calculus the equation becomes 
polynomial and leads to divisibility conditions of certain elementary polynomials. These 
can be solved for infinitely many orders at once. The necessary and sufficient equations 
for the ratio of eigenvalues are obtained directly without having to specify the nonlinear 
part explicitly. 

Assume that K and S are quadratic. A quadratic differential monomial is transformed 
into a symmetric polynomial in two symbols as follows 



The expression is symmetrized and divided by the number of symbol-permutations in 
order to ensure that 

ViVj = VjVi 

This procedure turns the operation of differentiation into ordinary multiplication 

DviVj = Vi+iVj + ViVj+i 

~ 2 

like the action of the Frechet derivative on a linear term 

Dv.VjVk = Vi+kVj + ViVj+k 

2 

The symmetry condition for quadratic polynomials {K,S € C[^i,^2]) reads 

Gn[ai,a2]S = Gm[bi,b2]K 

with the G-functions 

Gn[ai, a2](6, 6) =' ai(6 + 6)" - a2{^^ + i^) 

If Gm\bi,h'2\K is divisible by G„ [01,02] we have a symmetric polynomial expression for 
S which can be transformed back. Because the ^i-degree of K (the maximal number of 
x-derivatives of u in the terms of K) is smaller than n, the function Gn[ai,a2\ cannot 
divide K. Therefore G„ [01,02] should have a common factor with Gm[&i,^2] for there 
being a symmetry. Suppose we can find oi, 02, 61, 62 € C such that (F, L, T G C[^i, ^2]) 

G„ [01,02] = FL 
Gm[hiM]=FT. 

Then the Lie bracket (|2.l| ) vanishes if one takes K = LM and S = MT. One is free to 
choose M € C[^i,^2] as long as the ^i-degree of K remains smaller than n. If this sounds 
a bit too easy, the reader should note that this determines the system and its symmetry 
at the same time, and does not say anything about hierarchies of symmetries. 
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3 The use of resultants 

If the resultant of two polynomials vanishes, then their greatest common divisor has pos- 
itive degree. 

Lemma 1. Second order B-systems have symmetries at all orders. The ratio of eigenval- 
ues (and quadratic part) of the symmetries are fixed. 

Proof. Take ai ^ 02 again and r 7^ — 1. The G-function is 

G2[ai, 02] = ^ (6 - r^2)ir^i - ^2) with (02 - ai)r^ - 2air - (02 - ai) = 0, 

and the resultant of (^1 — r^2) with Gm[bi, 62] vanishes when 

61 _ 1 + r'" 
h ~ (1 + r)™' 

With this ratio (r^i — ^2) is a factor as well because the G-function is symmetric in ^1, ^2 
(the fraction ^j^^j^ is invariant under r — > i). ■ 

This implies that to cover systems with finitely many symmetries the degree of the common 
factor of the G-functioiis should be higher than 2. Degree 3 is not enough because a third 
degree symmetric polynomial always contains the factor + ^2- Their corresponding 
systems are always in a hierarchy of first, second or third order. 
We look at factors of degree 4. 

Lemma 2. The function Gn[l + r", (1 + r)"'](^i, ^2) has a factor of the form 

(6 - r-6)(r-ei - 6)(6 - s6)(s6 - 6) 

whenever Un =^ + r", (1 + r)"](l, s) = 0. 

Proof. The condition C/„ = is expressing the fact that the ratio of eigenvalues of the 
G-function containing root r equals the ratio of eigenvalues of the G-function containing 
root s. ■ 



In the following we disregard the trivial factors of C/„ which are (r — s){rs — 1) for all m 
and (r + l)(s + 1) when m is odd. 

Lemma 3. Take n > 3. To obtain all eigenvalues of n^^ order B-systems with a symmetry 
on order m one calculates the resultant of Un and Um with respect to s and applies the 
map r (\Xr)"^ roots. 

Proof. If the resultant of Un and Um vanishes for some r € C then by the previous Lemma 
Gn[l + r", (1 + r)"](6,6) and G,„[l + r™, (1 + r)'"](6,6) have a common fourth order 
factor. This implies that the n^^ order ,B-system with eigenvalues a\ = l+r", 02 = (1+r)" 
and quadratic part G„ [ai , 02] divided by this fourth order factor has a symmetry on order 
m. ■ 
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Example of Bakirov: The resultant of U4 and Uq with respect to s contains the factor 

/(r) = 2r^ + lOr^ + 15r^ + lOr + 2. 

We have that 

1 + r modj(r) = 5 

and 

(1 + 0^ mod /(r) . !l^!:!±|±M. 

Their ratio is 5, the ratio of the eigenvalues of the Bakirov system. As expected ^4(5, 1](1, r) 
is proportional to /(r). 

4 The use of p-adic numbers 

The use of p-adic methods in integrability theory was initiated in . For an introduction 
in p-adic number theory, see §. In this section we give a more expanded proof of the fact 
that the Bakirov system contains exactly one symmetry. 

The p-adic field is notated by Zp where p is some prime number. Its elements are 
represented by series of the form X^„>o OnP" with coefficients o„ G Z/p. The p-adic 
expansion of an positive integer is just its base p representation. For rational numbers we 
can get an infinite sequence. Examples: in Z5 we have 

57 = 2 • 5° -Fl • 5I 2 • 5^ 

| = l + T^ = 2-50 + Ei=il-5^ 

An element is invertible (in Z^ ) if it is nonzero modulo p, that is: oq 7^ 0. 
4.1 Hensels Lemma 

The following lemma gives a method to check whether a polynomial has a root in Z^ . 
Lemma 4 (Hensel). A polynomial 

m 

f{r) = Cjr* with Ci £ Zp 
i=0 

has a root in Z^ if there exists an a\ ^'Ljp such that 

• /(ai) = mod p 

• ^(ai) ^ mod p. 

Proof. It is possible to construct a sequence with 

On mod p"~^an-i 
f{an) = mod 
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Calculate [3 ^"L/p such that 



= /(a„+i) = /(a„ + = f{an) + -r{ai)l3p'' mod p 

ar 

By the induction hypotheses there exists a 7 G Z/p such that 

/(a„) = 7P" mod 

Substituting this and dividing by p^ gives an equation that can be solved in Z/p: 

Since the first step of the induction is part of the hypotheses, this concludes the proof. 
For example the square roots of 2 are in Z7. Take 

f(r) = r^-2 

Then 



/(3) = mod 7, /(4) = mod 7 
f (3) = 6 mod 7, f (4) = 1 mod 7 

So Hensels Lemma can be applied. The number 3 is lifted as follows. Modulo 7^ we have 



/(3) = 1-7 

so 7 = 1. The inverse of 6 in Z7 is 6. Then /? = —1 -6 = 1. Indeed 

/(3 + 1 • 7) = 2 • 7^ 

One step further gives 

/(3 + 1 •7 + 2-72) = 6-7^ +4- 7^ 
This shows that the method of Hensel is constructive. 

4.2 The method of Skolem 

Skolems method allows us to conclude that there exist only a finite number of symmetries. 
At first sight it looks a bit technical, but it is extremely powerful in our context. The 
method is based on the fact that if some equation does not have a solution in some p-adic 
field then it can not have a solution in C. Moreover the method reduces the number of 
orders that need to be checked to a finite number. 

If Xi G Zp then by Fermats little theorem there exists ayi & Zp such that x^~^ = l+Uip. 
Let 

1=1 

For instance, Un, as defined in lemma has the form of u'^ with Cj = (—1)* and j = 4. 
Lemma 5 (Skolem). If ^ mod p then Vr, u^_^_j.(^p_i^ / 0. 
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Proof. 

j 

ul+r{p-i) = Cixf (1 + VipY = ul mod p^O. 

i=l 

Therefore itself is 7^ 0. 

Lemma 6 (Skolem). //n^ = and uj^^O mod p then Vr > we have u^^^(^_i^ ^ 0. 
Proof. Assume tx^^^^p.^) = 



= J2 ^i^H^ + ViPf = E ( P'^'k 
1=1 t=i ^ 



use 







r — 


n 






. t- 





and divide by pr to get 



t=2 



This contradicts the second assumption since always contains a factor p. To see this 
write t = p°'s with p J( s. Then s is invertible and 

„t— 1 1 

The power of p is bigger than 1 for when a = we know s > 2 and when a 7^ we have 
s > 1 and > a + 2 (because p> 2). Hence we conclude u^^^^_i^ 7^ 0. ■ 

With the lemmas of Skolem one has to search a prime number p such that the Xi are in 
the field , and check the conditions for finitely many orders (p-2). The computations 
one has to do are all modulo p 01 p^. 

5 The Bakirov system 

Here is how to use these lemmas for the Bakirov system. We let p increase and look for 
p-adic roots of the resultant 2r^ + lOr^ + 15r^ + lOr + 2. The first prime such that all 
conditions are satisfied is 181. In Z/181 we find /(66) = /(139) = 0. These numbers can 
be lifted to elements of Z^g^ . Modulo they are 

r = 66 + 13p, s = 139 + 29p 

The function Um{r, s) is has the form with Cj = (— 1)', j = 4 and 

xi = 1 + s = 140 + 29p mod 
0:2 = 1 + r = 67 + 13p mod 

X3 = r(l + s) = 9 + 165p mod p^ 
X4 = s{l + r) = 82 mod p^ 
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For < m < 180 we have u^{r,s) = mod p only when m £ {0,1,4,6}. Applying 

x'r -1 . 
Xi — gives 

Hi = 40 mod p 
1/2 = 33 mod p 
7/3 = 140 mod p 
Hi = 46 mod p 

For m e {0, 1, 4, 6} the function 

33 • 66*" + 46 • 82™ - 40 • 140™ - 140 ■ 9™ 

is nonzero modulo p. Both Skolems lemmas can be applied and it is shown that there is 
no non-trivial symmetry but at order 6. 

6 The counter example to Fokas' conjecture 

Theorem 1. There exists a 2-component equation with exactly two non-trivial symmetries. 

Proof. The resultant of Uj and Uu has the following factor in common with the resultant 
of Uj and U29 

{r^ - r - l)(r^ + - l)(r'^ + 3r^ + 5r^ + 5r^ + 5r^ + 3r + 1) 

In Z/101 the first factor has solution 20 and the third solution 52. These can be lifted 
and both Skolems lemmas can be applied. In this way it is proven that the set {r, i, s, i} 
that corresponds to the modulo 101 set S = {20, 96, 52, 68} is no solution set of a Gm- 
function when m ^ {0, 1, 7, 11, 29}. If one of the other sets corresponding to {40, 48, 42, 89} 
or {32, 60, 63, 93} is a solution set of the function Gm for some m then their minimum 
polynomials divides the resultant of Uj with Um- That means that 5 is a solution set of 
Gm as well, hence m equals 0, 1, 7, 11 or 29. ■ 

We will compute the three equations and their symmetries explicitly. Each root of r^+r^— 1 
is mapped to a different eigenvalue. We take C[r]/(r^ + — 1) as our coefficient field. 
The eigenvalues of the systems will be 

1 + = 2r^ and (1 + rf = I6r'^ + 28r + 21 

Their quadratic part will be 

G7[2r^l6r2 + 28r + 21](6,6) 

divided by 

2(6 - -{r + + (1 - r - r^)^^ + ^2) 

i.e. 

!(2r2 + 4r + 3)(ei + + (2 - r)66 + Cl) 

Our examples, written more compactly than in [^], look like 

ut = 2r^UY + 7(2r2 + 4r + 3){v3Vo + (3 - r)v2Vi) 
vt = (16r2 + 28r + 21)?J7 



+ - 1 = 0. 
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The symmetries can be calculated in the same way, leading to 

ut= (-3r2 + r + 2)mii 

+11 ((14r2 + 24r + 18)t;7fo + {35r^ + 57r + 42)veVi 
+ (48r2 + 70r + 49)t;57;2 + (Slr^ + 65r + 42)v4V3) 

vt = (151r2 + 265r + 200)vii 

and 

ut = (-40r2 + 9r + 17)^29 

+29 (30(1081r2 + 1897r + 1432)t;25t;o 
+(311920r2 + 547311r + 413143)^24^^1 
+(706832r + 533441 + 403277r2)t;23V2 
+(449543r2 + 782050r + 589257)^;22^^3 
+(537572r + 402545 + 317304r2)/2f2it;4 
+(1026233r2 + 1635821r + 1205570)t;2oW5 
+(1101516r + 779787 + 787277r2)/2ui9U6 
+(2656229r + 1710194 + 2393075r2)wi8?;7 
+(3831912r2 + 3208669r + 1731205)?;i7?;8 
+(6105788r2 + 4007995r + 1678107)t;i6W9 
+(4807604r + 1421555 + 8899703r2)i;i5Uio 
+(5263833r + 11440843r2 + 915604)?;i4Wii 
+ 3(1793035r + 155000 + 4312473r2)?;i3Wi2) 

Vt = (3761840r2 + 6601569r + 4983377)^;29 

7 More symmetries 

We present the results of large computer calculations we did in MAPLE |Q] . We calculated 
the resultant of Un and Um for 4 < n < 10 and n + l<m<n + 150. To obtain the 
systems with finitely many symmetries only, one has to filter out the integrable systems. 
How to find all integrable n**^ order systems for fixed n if the quadratic part is Vq is 
described in This method we have extended to cover the case of an arbitrary quadratic 
part. This will be the subject of another paper. One finds common factors on orders 
m = k{n — 1), fen, {2k + l)n for all A; € N. 

To give an indication of the size of the expressions. The resultant of Uiq and f/ieo has 
degree 556. The coefficients of r" with 244 < n < 312 have 207 (decimal) digets. The 
number of n**^ order systems we have been calculating is 



n 


4 


5 


6 


7 


8 


9 


10 


4-10 


# 


2745 


2701 


5679 


5644 


8740 


8839 


11952 


46300 



In the pictures on the next pages the positions of the roots of these resultants in the 
complex plane are plotted. As a fundamental domain the upper half unit circle is chosen. 
The full pictures are invariant under r — > ^ and r — > f. 




Figure 1. Numerical approximations of the roots of tiie G-functions corresponding to almost 
integrable systems for order 4, 5 and 6 
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n=8 





Figure 2. Numerical approximations of the roots of tlie G-functions corresponding to almost 
integrable systems for order 7, 8 and 9 
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n=10 




Figure 3. Numerical approximations of the roots of the G-functions corresponding to almost 
integrable systems for order 10 



All these systems have at least one nontrivial symmetry. To answer the question how 
many symmetries there exactly are we implemented the method of Skolem. We made the 
following refinements. 

• Most of the resultants we have calculated are irreducible. By the argument in the 
proof of theorem ^ it suffices to prove the statement for one particular set of roots. 

• Sometimes it is much more efficient to use two pairs of roots. The argument goes as 
follows. The resultant of and Uig contains the factor 

f[r) = _^ 4^11 _^ iQ^io _^ iQ^9 _^ 28r^ + 34r^ 

+37r^ + 34r5 + 28r^ + I9r^ + lOr^ + 4r + 1 

which is irreducible over Q and splits into linear factors over Z^qq. The numbers 
(264,407) are a solution for Um{r,s) when m £ {0,1,5,19,256,414}. The numbers 
(267, 300) are a solution for Um{r, s) when m G {0, 1, 5, 19, 162, 254}. By using both 
pairs we can apply lemma |5| for all < m < 508 but {0, 1, 5, 19}, for which we can 
use the lemma ^. The computer could not find any prime such that the normal 
procedure works, it has been busy for days to check all primes p < 8147. 

With these improvements we have been able to prove that all these system have exactly one 
non trivial symmetry, with the exeption of the seventh order systems with two symmetries 
at order 11 and 29. 

The following MAPLE output can be used to verify the above statement for n = 7, 29 < 
m < 37. 



prf29: = [101, {20, 52}] , [97, {4, 32}] 
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prf30:=[2531, {75, 871}]: 
prf31:=[1021, {16, 42>] : 
prf32:=[877, {226, 214}]: 
prf33:=[601, {23, 409}] : 
prf34:=[2857, {2457, 716}, {742, 391}]: 
prf35:=[661, {401, 330}, {122, 245}]: 
prf36:=[179, {17, 76}] : 
prf37:=[233, {30, 56}, {20, 84}]: 

The sequence prf . m contains the proofs for different factors of the resultant of Uj and 
Um- Each proo/ consists of a prime number p and one or two sets of modufo p solutions 
such that all conditions of Skolem are satisfied. The exceptions, where the resultant has 
two factors, are 

(n, m) = (4, 24), (4, 28), (6, 42), (7, 8), (7, 49), (8, 56), (10, 70) 

Three factors appear at n = 7, m = 11 and four at n = 7, m = 29. 

8 Conclusions 

We have shown that the existence of one or two generalized symmetries of an evolution 
equation does not necessarily imply integrability. We hope that this illustrates the use 
of p-adic and resultant methods to this field and that these methods will be more widely 
applied. With these results in mind this puts a burden of proof on anyone claiming 
integrability (with respect to generalized symmetries). We mention the successful use of 
number theoretic methods, especially the Lech-Mahler theorem, in this respect, cf. |T^. 
These methods are not restricted to the special kind of systems we study here, but they 
are applicable to any polynomial system, in principle. 
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